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The contribution of the -R-parity violating supersymmetric model to the fermion electric dipole 
moment at the two-loop level is analyzed. We show that in general, the Barr-Zee type contribution 
to the fermion electric dipole moment with the exchange of W and Z bosons is not small compared 
to the currently known photon exchange one with R-parity violating interactions. We will then give 
new upper bounds on the imaginary parts of .R-parity violating couplings from the experimental 
data of the electric dipole moments of the electron and of the neutron. 
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I. INTRODUCTION 

The standard model (SM) of particle physics, although 
being very successful in interpreting many experimental 
data, has difficulty in explaining some phenomena such 
as the matter abundance of our Universe. The SM has 
therefore to be extended. 

There are currently many approaches to search for new 
physics (NP) beyond the SM. Among many others, the 
measurement of the electric dipole moment (EDM) is of 
particular interest. The reasons are the following. The 
EDM is an observable sensitive to the violation of the 
parity and the time-reversal symmetry (or equivalently 
the CP). The contribution from the SM is in general 
very small and this fact makes the EDM a sensitive 
observable to NP with large CP violation. The experi- 
mental data available are very accurate for a variety of 
systems such as the neutron (d n < 2.9 x 10 _26 ecm) 0, 
the 205 T1 (d T i < 9 x l(T 25 e cm) 3] and the 199 Hg atoms 
(d Hg < 3.1 x l(T 29 ecm) 4], the YbF molecule (which 
ives bound on the electron EDM: d e < 1.05 x 10~ 27 e cm) 
the muon @, etc. The EDM is therefore a very good 
probe of NP. New generation of experiment using storage 
rings is also in preparation, aiming at the measurements 
of the EDMs of muon, proton and deuteron Q. 

On the theoretical side, the minimal supersymmetric 
standard model (MSSM) is known to be the leading can- 
didate of NP Q . A general supersymmetric extension of 
the SM allows baryon number or lepton number violat- 
ing interactions, so we must impose the conservation of 
R-parity (R = (_l)3B-£+2s) to forbid them _ This as _ 

sumption is, however, completely ad hoc, so the .R-parity 
violating (RPV) interactions have to be investigated phe- 
nomenologically. Until now, many RPV interactions were 
constrained by high energy experiments, low energy pre- 
cision tests, and cosmological observations @. Thanks 
to many efforts in EDM experiments, many phenomeno- 
logical analyses of models of multi-Higg s fl0i - [l3| , super- 
symmetry with [T^ - [r8| and without [l9l427l | .R-paritv in- 



variance were done, and many CP phases have been con- 
strained so far. 

Previous analyses showed that RPV interactions con- 
tribute to the fermion EDM starting from the two-loop 
level [2(j, and that the Barr-Zee t ype t wo-loop level di- 
agrams give the leading effect [2l|, (25|. It was argued 
that the fermion EDM receives the largest part from the 
photon exchange Barr-Zee diagram, and that the other 
Barr-Zee type diagrams with W and Z bosons are sub- 
leading, on the base of the analogy between the .R-p arity 
violation and models with charged Higgs boson [12|, [21j . 
We must be careful however to the fact that these two 
processes cannot be described similarly, since the SU (2) l 
gauge structure (chirality structure) of the Yukawa inter- 
actions with charged Higgs boson differs from that of the 
.R-parity violation. As the chirality structure inside the 
fermion loop of the Barr-Zee type diagram is different, 
the naive analogy no longer works. Moreover, the RPV 
Barr-Zee type diagram with W boson exchange changes 
the flavors of the fields in the loop, and consequently 
many additional combinations of RPV couplings with 
new flavor structure may be constrained from the the 
EDM experimental data. 

The main purpose of this paper is then to evaluate 
the two-loop level Barr-Zee type diagram with W and 
Z boson exchange within .R-parity violation, and to up- 
date the constraints on the RPV interactions provided by 
the current experimental data of the EDMs of the YbF 
molecule, the neutron and the 199 Hg atom. 

This paper is organized as follows. We first present the 
RPV interaction in the next section. We then present and 
formulate in Section [TTTI the two- loop level Barr-Zee type 
diagram with Z boson exchange. In Section lTVl we derive 
the formula for the RPV Barr-Zee type contribution with 
W boson exchange. This will be done in two steps, first 
by constructing the inner fermion loop. We then attach 
this loop to the external fermion line. In Section [Vj we 
analyze the limits provided by the current experimental 
data of the YbF molecule, the neutron and the 199 Hg 
atom. We finally summarize our work in the last section. 
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II. RPV CONTRIBUTION 

The superpotential of the RPV interactions relevant in 
this discussion can be written as follows: 



1 



+ \K jk {U c )i{D c ) j {D c ) k , 



(1) 



with k = 1,2,3 indicating the generation, a,b = 1,2 
the SU(2)l indices. The SU(3) C indices have been omit- 
ted. L and E c denote the lepton doublet and singlet 
left-chiral superfields. Q, U c and D c denote respectively 
the quark doublet, up quark singlet and down quark sin- 
glet left-chiral superfields. The bilinear term has been 
omitted in our discussion. We also neglected the soft 
breaking terms in the RPV sector. Also baryon number 
violating RPV interactions (A^ fc ) will be omitted from 
now, to avoid rapid proton decay. This RPV superpoten- 
tial gives the following lepton number violating Yukawa 
interactions: 
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The projection of the chirality is given by Pl = 5(1—75) 
[and we also define Pr = ~(1 + 75) for later use]. 



III. BARR-ZEE TYPE DIAGRAM WITH Z 
BOSON EXCHANGE 

The RPV interactions contribute to the fermion EDM 
starting from the two-loop level. In this section, we will 



give the formula for the Barr-Zee type diagram with Z 
boson exchange (see Fig. [TJ. The computation of this 




FIG. 1. An example of Barr-Zee type diagram with Z boson 
exchange generated with RPV interactions. 

diagram is very similar to that of the Barr-Zee type dia- 
gram with photon exchange [25| . 

The first step of the evaluation of the Barr-Zee type 
two-loop diagram is to construct the gauge invariant ef- 
fective v-fZ vertex with the fermion loop diagrams, as 
shown in Fig. [5J This method is based on the analyses of 
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FIG. 2. One-loop effective DyZ vertex generated with RPV 
interactions. 



Refs. [TJ, [13|, |21| . We will then attach this gauge invari- 
ant effective vertex to the external fermion line to obtain 
the EDM operator. The amplitude of the inner fermion 
loop is given as 



iMo 1 z = -^i 33 n c Q f e 2 a f e* (qi)el(q 2 ) 
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where the last line was given by taking the leading order 
contribution of the expansion in q\ . This approximation 
is justified since the EDM is the first order coefficient 
of the multipolar expansion. Here % and j are the fla- 
vor indices of the incident sneutrino and loop fermion /, 
respectively, and A is the RPV coupling, where A = A 



for inner loop charged lepton, and A = A' in the case of 
quark. For quark loops, the number of color is n c = 3 
(for lepton loop, n c — 1). The weak coupling is given by 
ctf (f = I, d), where ai = |(3tan6>w — cotOw) ~ —0.065 
for the coupling of the Z boson with charged leptons, 
and otd = J2 tan^w — \ cot 9w ~ —0.42 for the coupling 
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with quarks. We must note that a/ is the vector coupling 
of the Z boson with fermions, and that the axial vector 
coupling does not contribute. 

The next step is to put the above effective v^Z vertex 
into the external fermion line to form the EDM operator. 
This manipulation is again very similar to that for the 
photon or gluon exchange Barr-Zee type diagram [25| . 
and is given independently of the gauge as follows: 
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where /(a, b, c) is defined as follows: 
I{a,b,c) = 
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In the last approximation of Eq. we have used 

to~. >> m 2 z » rrvi.. Note that the fermions fj and 
Fk are either down-type quark or charged lepton. The 
RPV coupling A = A for external charged lepton, and 
A = A' for the case of quark. 

The EDM cLf is defined as follows 



-EDM 
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where F M „ is the electromagnetic field strength. The 
fermion EDM with Z boson exchange generated by RPV 
interactions is then given by 
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Let us compare the Z boson exchange contribution with 
the photon exchange process. The photon exchange con- 
tribution is given by [25| 
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We remark that dp and d Fk of Eqs. (JT)) and (HJ have the 
same sign. The loop factor (the logarithmic factor) has 

ml 

the same order of magnitude, with In — ^ being around 
( m f\ 

one half of — 2 + In — j 2 - for sneutrino masses of or- 
der TeV. For lepton EDM with lepton loop, lepton EDM 



with quark loop and quark EDM with lepton loop, the 
Z boson exchange contribution is small compared to the 
photon exchange one, since the weak charge of the lepton 
is small (a; = —0.065). For quark EDM with quark loop 
however, both diagrams have the same order of magni- 
tude. This means that an important enhancement of the 
RPV contribution will occur. The detailed analysis will 
be done in Section IVl 



IV. BARR-ZEE TYPE DIAGRAM WITH W 
BOSON EXCHANGE 

The computation of the Barr-Zee type diagram with W 
boson exchange is similar to that of the Barr-Zee type 
diagram with charged Higgs exchange, as was done in 
Ref. [IH ■ Here we should give the detail of the derivation. 



A. Inner fermion loop 
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FIG. 3. One-loop effective Sl^/W vertex generated with RPV 
interactions. Note that for the lepton loop, diagram (a) does 
not contribute, since the neutrino has no charge. Diagram 
(d) is the Nambu-Goldstone boson (cj>) exchange contribution, 
which disappears in the non-linear R$ gauge. 

As for the previously seen Z boson exchange contri- 
bution, we first derive the gauge invariant expression for 
the one-loop effective e^W vertex, generated by RPV 
interactions. The contributing diagrams are shown in 
Fig. [3] In our calculation, we have chosen the non-linear 
i?5 gauge [11] which is given by the following gauge fixing 
(W boson) 



J -GF 



(9) 



In this gauge, the calculation becomes easy (interaction 
between Nambu-Goldstone boson (</>) and gauge boson 
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cancels) . The derivation of the one- loop ej/yW amplitude 
goes in a manner very similar to that of the fermion-loop 
contribution to the charged Higgs boson decay into W 
boson and photon (29[. 

The amplitude of the one- loop effective e^jW vertex 



generated with RPV interactions is given by 
where 
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where i, j and a denote respectively the flavor of the incident selectron, down type and up type quarks of the loop. 
The convention for A and A are the same as for the previous case with Barr-Zee type diagram with Z boson exchange. 
The labels /" and f d denote respectively the up and down type quarks for quark loop, the neutrino and charged 
lepton for lepton loop. For the quark loop, the number of colors is n c = 3, and n c = 1 for lepton. Moreover, we 
have Q u = and Qd = — 1 for the lepton loop, and Q u = | and Qd = — \ for the quark loop. Here V a j is the CKM 
matrix element for the quark loop contribution. For the case of lepton loop, V a j is a simple unit matrix, as we do 
not consider the flavor mixing in the lepton sector. It should be noted that the diagram (d) in Fig. [3] vanishes in the 
non-linear i?j gauge. The second and third lines of Eq. (fT5")) become 
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312 



independently of the constant f, where q\ 2 — q\ + q%. 

The calculation of the one-loop integral can be performed by using the Passarino-Veltman one-loop tensor [3(| (see 
Appendix [X] for detailed derivation) . By taking the leading gauge invariant contribution first order in q\ , we obtain 
the following amplitude: 
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Diagrams (a), (b) and (c) are each divergent, but the di- B. Second loop 

vergence of the total contribution cancels out, as is shown 
in Eq. ([Mil . 

The Barr-Zee type EDM with W exchange can be con- 
structed by attaching the effective one-loop level e^W 
vertex into the second loop as shown in Fig. 0] The total 
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RPV Barr-Zee type contribution with W boson exchange is given by 
I 
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where we have kept only terms contributing to the EDM, 
and we have used the formula of Eq. ([5]). We remark that 
the second term in the curly brackets 

-z [( qi ■ k')g^ - k'^ql + ie^( qi ) a k' ] 

in the z integration vanishes. This cancellation is rem- 
iniscent of that we encountered for the RPV Barr-Zee 
type EDM with exchange of photon and gluon [25| . 

We thus obtain the next formula for the EDM with W 
boson exchange: 
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FIG. 4. Barr-Zee type contribution to the fermion EDM with 
W boson exchange generated from RPV interactions. The di- 
agram on the left side involves the complex conjugated com- 
bination of RPV couplings. F u and F d denote respectively 
the up and down type quarks for quark EDM, the neutrino 
and charged lepton for lepton EDM. 
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We treat first the case where there is no top quark in 
the loop. In this case, the masses of quarks and leptons 
can be neglected since rn 2 , u ,m 2 d << m 2 v ,ra\ L .. The 

EDM is then given by 
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where we have assumed that °m\ Li » m 
n c (Q u + Qd) = +1 for inner quark loop and s/ = —1 for 
lepton loop. We have verified numerically that the above 



approximation works well when top quark is absent in 
the loop. 

For the case where top quark is present in the loop, 
the mass m 2 „ cannot be neglected. The integral of Eq. 
(|17[) can be performed analytically, and we obtain 
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where z\y = m^/m^ and zg = m~ L ./mf, and L'i2(z) is 
the dilogarithm function (for the calculation of the dilog- 
arithm function, we have used the computational code of 
Ref. (HI). The integral with top quark in the loop is of 
course smaller than the case without it. For mg £j = 1 
TeV, the loop integral with top quark becomes about the 
30% of that without top quark. 



V. ANALYSIS 
A. Constraints on Im(\ijjX* kk ) 

Let us first discuss the case where no change of genera- 
tion occurs in the inner fermion loop and external fermion 
line (i.e. a — j and b = k). This contribution interferes 
with the Barr-Zee type process with photon orgluon ex- 
change, previously investigated in Refs. |^ |20h23 . |24| . |25| . 
This part is thus an extension of these analyses. In this 
analysis, we assume that all sparticle masses are equal to 
1 TeV [33. 

The first case to treat is the lepton EDM with inner 
lepton loop, where the combinations of RPV couplings 
Im(A3nA3 22 ) and Im^iiA^) are involved. This con- 
tribution can be constrained by the current experimental 
data of the electron EDM given by the YbF molecule ex- 
periment 

|4| < 1-05 x KT 27 ecm . (20) 

The Barr-Zee type lepton EDM with lepton loop receives 
negligible contribution from the Z boson exchange dia- 
gram, since the charged leptons have small coupling with 
the Z boson (a e sa —0.065). The W boson exchange ef- 
fect is also small compared to the photon exchange one 
by about one order of magnitude. We must also note that 
the sign of the W boson exchange Barr-Zee type EDM 
is the same as the photon exchange EDM. The photon 
exchange contribution thus dominates, as was claimed in 
previous works. The upper limit to the combinations of 
RPV couplings is shown in Table Q] To our best knowl- 
edge, the experimental data of the electron EDM given 
by the YbF molecule experiment give the tightest limits. 

The second case to discuss is the lepton EDM with 
inner quark loop. In this case also, the Barr-Zee type 
diagrams with Z boson exchange is small, due to the 
small weak coupling. Numerically, it is smaller than 5%, 
with the same sign as the photon exchange EDM. The 
W boson exchange contribution is however not negligi- 
ble. This is because the photon exchange EDM receives 
a suppression by the fractional charge. For the bot- 
tom quark loop contribution, the ratio is /dj ss 0.25. 
Moreover dj and d^ have opposite sign. This means 
that the total EDM contribution becomes significantly 
small. By using the experimental data of Eq. ([20)) , it 
is possible to constrain the combinations of RPV cou- 
plings Im(AaiASi) (i = 2,3), Im(A m Ag 2 ) (i = 1,3) 
and Im(AiuAi3 3 ) (« = 1>2). For these RPV bilinears, 



TABLE I. Upper bounds to the RPV couplings given by the 
current EDM experimental data of the electron and neutron, 
with ma i = mg Li = 1 TeV. For comparison, we have also 
written the limits provided only by the photon exchange Barr- 
Zee type contribution. Limits given by other experiments 
[9I |26|] are also shown, i = 1,2, 3. 

RPV couplings This work d 1 only Other 
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there already exist stronger constraints, given by the ex- 
perimental data of the 199 Hg atom [1 via P, CP-odd 
electron- nucleon interaction [22|, HH, |34| . It is not possi- 
ble to give new upper limits on the corresponding RPV 
interactions with Barr-Zee type process. 

The third case to consider is the quark EDM with lep- 
ton inner loop. The Z boson exchange contribution has 
a small contribution for the same reason as the previous 
case (numerically, less than 5%). The W boson exchange 
EDM is however sizable, with about 38% of the photon 
exchange EDM for r loop contribution, and 24% for /1 
loop contribution. The quark EDM with lepton loop has 
the same sign as the photon exchange Barr-Zee diagram, 
so d w and d 1 interfere constructively. We obtain thus 
tighter limits for the RPV couplings. By using the cur- 
rent experimental data of the neutron EDM [2] 

\d n \ < 2.9 x 10" 26 ecm . (21) 

and the relation between the neutron and quark EDMs 
calculated with the QCD sum rules [35|, [36| 

d n = 0.47d d - 0.124 + e(0.354 + 0.174) , (22) 

the combinations of RPV couplings Im(Ai22A^* 1 ) (i — 
1,3) and In^A^A-^) (i — 1, 2) can be constrained. The 
result is given in Table |TJ In this case we could give a new 
upper limits on Im^A^A^) (i = 1,2). 

The final case to consider is the quark EDM with inner 
quark loop. In this case both Z and W boson exchange 
contributions are sizable compared with the photon ex- 
change one. The first reason is that the photon exchange 
EDM receives a double suppression from the fractional 
charge of the down type quark. For the Z boson exchange 
process, the suppression due to the small weak coupling 
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between Z and charged lepton is absent so that d z and 
d 1 have similar size. The sign between them is the same 
and the interference is constructive. For the W boson, 
the contribution is around 35% of the photon exchange 
EDM. We must note that the Barr-Zee type diagram with 
W boson and inner quark loop act destructively against 
the rest. The RPV interactions contributing to the quark 
EDM with quark loop can be constrained by the neutron 
EDM experimental data. These RPV interactions are, 
however, already strongly constrained by the experimen- 
tal data of the 199 Hg atom Bvia chomo-EDM and P, 
CP-odd 4-quark interactions [33L l35l H3] ■ It is not possi- 
ble to give new upper limits on the corresponding RPV 
interactions with Barr-Zee type process. 

B. Constraints on Im(V a j A* 6fc ) 

The Barr-Zee type diagram with W boson exchange 
provides also possibilities to constrain new combinations 
of RPV couplings through the generation change of the 
CKM matrix. By combining Eqs. {TSJ and jJSJ with the 
EDM experimental data (|20[) and (|2"Tj) (with the relation 
(|22[)). we obtain upper limits on the imaginary parts of 
RPV couplings shown in Table [TTJ 

We see that the experimental data of the electron 
EDM give new upper limits on lln^A^nA^)! and 
|Im(A3 11 A32 3 )|. These upper bounds could be obtained 
thanks to the enhancement due to the large bottom quark 
mass in the inner loop. Other limits on RPV couplings 
given by this analysis are weaker than those obtained 
by other experiments, but we remark that many bounds 



on RPV couplings are not so different. The progress on 
neutron and atomic/molecular EDMs are very promis- 
ing, and many combinations of RPV couplings are in the 
reach of the next generation experiments. The Barr-Zee 
type EDM with W boson exchange thus provides a very 
wide accessibility to the CP violation of the RPV sector. 

VI. CONCLUSION 

In conclusion we have discussed additional contribu- 
tion to the fermion EDM in the RPV supersymmetric 
models. We have calculated the contribution of the two- 
loop level Barr-Zee type diagram with W and Z bo- 
son exchange. We have then found that these contri- 
bution are not negligible in many situations, contrary to 
the claim in Ref. [21(. In particular, the quark EDM 
with lepton inner loop was enhanced with the W bo- 
son exchange diagram and we have done a significant 
update of the upper limits to the RPV couplings as 
|Im(A*33A^ n )| < 2.5 x 1CT 2 (j = 1,2). We have also 
found that many additional RPV contributions gener- 
ated by the flavor change due to the W boson exchange 
exist, with many combinations of RPV couplings. In this 
case it was also possible to give a new constraint on RPV 
couplings as |Im(A* n A^ 23 )| < 1.6xl0" 2 (j = 2,3). More- 
over, many RPV couplings can potentially be constrained 
by near future EDM experiments, since the upper limits 
given by this analysis and those from other experiments 
have close value. The W and Z boson exchange contri- 
bution thus provides a very wide possibility to approach 
the CP violation of the RPV sector. 



Appendix A: Passarino-Veltman one-loop integral and the El^W amplitude 



For the calculation of the one-loop e-LilW process (Eqs. (fTTj) , (jT2J) and (113[1 ). it is convenient to use the formalism 
of Passarino-Veltman one- loop integral [30| • 

[fc 2 — 771q] [(fc + pi) 2 — m 2 ] ' ^ ^ 

d k [k 2 - m i}[{k+p l y~ m \] ' (A2) 

[fc 2 - ml] [(k +pi) 2 - m\] [(fc +P2) 2 - m 2 ] ^ 

[fc 2 - ml] [(fc + Pl ) 2 - m\] [(fc + p 2 ) 2 -m 2 ] 1 1 ' 

^ r,., ,,, 2ir „.../*L, 1 r,,....„ =51 • ( A5 ) 



Let us define the one 


-loop tensor as 






Bo(pl,ml,mf 


_ (2^f 




£^(p 2 ,m 2 ,™ 2 ; 


, _ (27T/X) £ 
17T 2 


Co(pf, (pi - 


\2 2 2 2 2' 
Pa) ,P 2 ,m ,m 1 ,m 2 


, _ (^r 

in 2 


C^pKpi- 


\2 2 2 2 2' 

P2j ,p 2 ,m ,m 1 ,m 2 


, _ (2^) £ 
iir 2 


C^(pl,(p x - 


\2 2 2 2 2 

Pa) ,P 2 ,m ,m 1 ,m 2 


, _ (2^) £ 
in 2 



where d = 4 — e is the space-time dimension shifted by e. Due to the Lorentz covariance, B^, C M and C^ v can be 
decomposed as follows 

B» = B lP » , (A6) 
C» = + CzpZ , (A7) 
C» v = C 00 g» v + Cutfrt + C n (rtpl + p»p\) + C 22 p^ 2 , (A8) 
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TABLE II. Upper bounds to the RPV couplings given by the current EDM experimental data for mg Lj = 1 TeV. Limits from 
other experiments [§. are also shown. 



RPV couplings 


Limits given by this analysis 


Other experiments 


lm(A 2 iiA22i j| 


2.5 (e EDM [5]) 


2.0 x 10~ ; ' 


|Im(A 311 A 3 2i)| 


2.5 (e EDM [5]) 


2.0 x 10" 5 


|Ini(A 2 iiA 2 3i)| 


210 (e EDM [5]) 


8.2 x 10~ 4 


|Im( A 311 A 331 ) | 


210 (e EDM [5]) 


8.2 x 10" 4 


it ( \ * \ ' ^ i 
|lm^A 2 nA2i2j| 


0.12 (e EDM [5]) 


3.3 x 10~ 3 


IT { \ * \ ' '\ 1 
l im (A311 A 312,)| 


0.12 (e EDM [5J) 


3.3 x 10" 3 


IT ( \ * \ ' A 
|lm ^211^232 ) \ 


2.1 (e EDM [5]) 


2.9 x 10" 2 


|Im( A 311 A 332 ) | 


2.1 (e EDM [5]) 


2.9 x 10" 2 


|Im(A 2 nA2i3)| 


0.18 (e EDM [5]) 


2.9 x 10~ 2 


|Iin( A 311 A 313 ) | 


0.18 (e EDM [5]) 
1.6 x 10" 2 (e EDM [5]) 
1.6 x 10" 2 (e EDM [5]) 


1.7 x 10" 2 


|Im(A 2 iiA223)l 


2.9 x 10" 2 


|Im(A 311 A323) 


1.7 x 10" 2 


Im(A 12 2A 1 2i)| 


6.6 (n EDM [2]) 


2.9 x 10~ 2 


|Im(A 32 2 A 32 i) | 


6.6 (n EDM [2]) 


2.9 x 10" 2 


|Im(Ai22 Ai 3 i) | 


180 (ra EDM [2]) 


0.70 


|Im(A 32 2A 331 ) 


180 (n EDM [2]) 


0.70 


Im(A 133 Ai2i)| 


0.39 (n EDM [2]) 


1.7 x HT* 


Im(A233A 2 2l)| 


0.39 (n EDM [2]) 


2.9 x 10~ 2 


Im(A 133 A 131 )| 


11 (n EDM [2]) 


0.42 


Im(A233A 2 3l)| 


11 (n EDM [2]) 


0.70 


|Im(A-nA- 12 )| (i = 1,2,3) 


7.0 (n EDM [2]) 


7.3 x 10 -4 


MA£iAj sa )| (J = X ' 2 ) 


130 (n EDM [2]) 


8.0 x 10" 2 


|Im(A 3 * 11 A 3 32)| 


130 (ra EDM [2]) 


1.4 x 10~ 2 


|I m ( A l*21 A 112) 


31 (n EDM [2]) 


4.0 x 10" 4 


|Im(A' fc21 A^2)| (fc = 2,3) 


31 (n EDM [2]) 


3.2 x 10~ 3 


|Im(A-2iA-2 2 )| (i = 1,2,3) 


7.1 (n EDM [2]) 


7.3 x 10" 4 


|Im(A-2iA- 32 )| (i = 1,2,3) 


550 (n EDM [2]) 


8.0 x 10~ 2 


|Im(Ai 31 Ai 12 )j (i = 1,2,3) 


800 (n EDM [2]) 


8.0 x 10" 2 


l lm (-^31 A i22)l {i = 1,2,3) 


190 (n EDM [2]) 
1.5 x 10 4 (n EDM [2]) 


8.0 x 10" 2 


|Im(Ai 31 A^ 32 )| (i = 1,2,3) 


7.3 x 10" 4 


IMXIjA^)! (j = 1,2,3) 


11 (n EDM [2]) 


3.2 x 10" 3 


|Im(A-HA- 23 )| {i = 1,2,3) 


0.93 (n EDM [2]) 


3.2 x 10" 3 


|Im(A£ 21 A£ 13 )| (i = 1,2,3) 


47 (n EDM [2]) 


3.2 x 10~ 3 


llmfA'Si A'oq)! (j = 1,2,3) 


4.0 (n EDM [2]) 


3.2 x 10" 3 


|Im( A^2l -^133 ) I 


0.53 (n EDM [2]) 


2.0 x 10" 4 


|Im(Afc*2i Afc 33 ) 1 (fc = 2,3) 


0.53 (ra EDM [2]) 


8.0 x 10 -2 


|Im(A 131 A 113 )| 


1.3 x 10 3 (n EDM [2J) 


1.7 x 10~ 5 


|Im(A' fc * 31 A^ 13 )| (fc = 2,3) 


1.3 x 10 3 (n EDM [2]) 


8.0 x 10 -2 


|Im(Ar 31 Aj 23 )| (i = 1,2,3) 


110 (n EDM [2]) 


8.0 x 10~ 2 


|Im( A ^ 31 A J 


14 (n EDM [2]) 


4.9 x 10 -3 


|Im(A^ 31 A' fe 3 3 )| (fc = 2,3) 


14 (n EDM [2]) 


2.0 



where the arguments of the loop functions were omitted. The loop functions Bq, B\, Co, Coo, Cn, C12 and C22 can 
be explicitly calculated, but they satisfy also many useful relations. Note also that Bq, B\ and Coo are divergent. By 
contracting C with external momenta Pi or p% , we obtain 



ml + ml] C , 

m 2 2 + ml] Co . (A9) 



C »P\ = g (P2 > «■&<)> m l) - 2 B o((pi - V2) 2 ,m\,m\) - - [pi - 
C i*P% = 2 B o(Pi,ml,ml) - -B ((pi - p 2 ) 2 , m\, mj) - -[p%- 



Here we have omitted the arguments for Co which are implicitly the same as those of Eq. (|A3[) . 
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Similarly, we obtain also additional relations by contracting C^ u with two external momenta. 



dC 00 + Gup\ + 2C\ 2 {pi -Pi) + C22P2 = B a{{P2 - Pi) 2 ! ,mf,m!) + TOqCo , 

Coo + CupI + C 12 {pi -p 2 ) = ^i((p 2 -pi) 2 ,m?,m,2) + ^B ((p 2 - pi) 2 ,m\,m\\) 



C12P1 + C 22 (pi -p 2 ) 



~2 [pi~ m i + m o] Ci > 

2-Bi(p2,mo,7n|) - ^-Bi((p 2 - pi) 2 , m 2 , m 2 .) 
1 



Cii(pi • p 2 ) + C12P2 = 



~2 bi- m i+ TO o] C2 , 
2-Bi(pi,mo,mi) + --^((Pa - pi) 2 , m 2 , m 2 ) 
+ 2 B o((P2 -pi) 2 ,mi,TO2) ~ o _TO 2 + TO o] Ci > 
Coo + C 12 (pi -p 2 ) + C22P2 = -^B 1 ((p 2 -pxf,m\,ml) - - [pi - m\ + to 2 ,] C 2 , 
where again arguments for C functions were omitted. Applying all these relations to Eqs. (fTTj) . we obtain 



(A10) 
(All) 

(A12) 

(A13) 
(A14) 



«M( a ) = — — T^- A 



™ C e*(<Zi)e*(g2) 



(4tt) 2 ^y^sin^' 

x { [((ft • «a)«r - <#tf)(2C ia + C 2 ) + ie^(9i) a fe)^C 2 ] (g 2 , ( Ql + q 2 ) 2 , g 2 , to 2 „, m 2 „, to 2 ,) 

+5i((gi+g 2 ) 2 ,m 2 ? ,TO 2 / ) 5 ^} . (A15) 

Note that the expression in parenthesis (q 2 , (qi + q 2 ) 2 ,q 2 ,m 2 fu ,m 2 f u,m 2 , d ) denote the common arguments for the loop 

> i 3 Jj 

functions C\ 2 and C 2 . For Eqs. (fTS"]) and (JT3J) , we obtain similarly 



iM(b) = 



iM {c) = 



2im f d Oip 2 V ■ 
(4n) z y/2sm6w 

x {[((91 • q2)9^ - g£tf)(2C ia - C 2 - Co) - ie^( qi ) a (q 2 ), 3 {C 2 + Co)] (qf, (gi + g 2 ) 2 , q%, m) f , m) f , m}u) 
B 1 {{q 1 +q 2 ) 2 ,m 2 f d,m 2 fr ) + B ((qi+q 2 ) 2 ,m 2 f d,m 2 fr )'] g»"} , (A16) 

Tr -'Ha J ~/ Va ' n ce ;( gl )e:( g2 )ff^i3 1 ((g 1 +g2) 2 ,m 2 ,,m 2 d ) . (A17) 
(47r) 2 y / 2sm9 w 1 3 



2im f d 

jj 



We should note that the arguments of the C loop functions for iM. (t>) differ from those for iM. ( a ) . 

We should now show that terms with B functions cancel with each other for iM^ + iM(h) + tM.t c y The loop 
functions Bq and B\ can be directly calculated using Feynman parameters as 



B (q 2 ,ml,ml) 
Bi(q 2 ,ml,mf) 



7 + ln(47r) 



- 7 + ln(47r) 



dx In 
■1 



-x(l — x)q 2 + (1 — x)mQ + xm\ 



dx x In 



-x(l — x)q 2 + (1 — i)mg + xm\ 

7 2 



(Ais) 

(A19) 



where 7 is the Euler constant and [i the mass scale peculiar in the dimensional regularization. By noting that 



B Q (q 2 1 ml,ml) + ^(g 2 , to 2 ,, m 2 ) = - 



7 + ln(47r) 



dx (1 — x) In 



-x(l — ir)g 2 + (1 — .t)toq + xm\ 

7 



7 + ln(47r) 



dx x In 



-x(l — x)q 2 + (1 — x)m\ + xm 2 



-5i(g 2 ,m 2 ,m 2 ) 



(A20) 
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the sum of the B functions for iAi( a ) + iM-Oo) + iM-(c) becomes 

Q u Bi{(Qi + q 2 ) 2 ,m 2 ^,m 2 fd ) + Qd + q 2 ) 2 , m 2 d , m 2 u) + B ((qi + q 2 ) 2 , m 2 d , mju) 

-Bi((qx + q 2 ) 2 ,m 2 u,m 2 d ) 
= (Qu - Qd - l)Bi((qi + q 2 ) 2 ,m 2 ^,m 2 fd ) 
= . 



(A21) 



We see that the divergence cancels. The evaluation of the remaining C functions can also be done by direct calculation 
of the integrals. 



1 n-\-Z 

dz dx 



Ci2(<?i, (qi + q2) 2 -,ql,m 2 f u,m 2 fu ,m 2 fd ) 

J J J 7" 



C2(<?i,(<7i +92) > q 2 ,m f u,m f u,m fd ) = dz 

^ Jo Jo 



xz 



o x(l - x)q{ + z(l - z)q% + 2xz(qi ■ q 2 ) - (1 - z)mi u - zm fd 

j J j 

(A22) 



1 l-l-Z 

dx 



Ci 2 (q 1 ,(qi+ q 2 ) , q 2 ,m fd ,m fd ,m f u) = j dz 

s~t f 2 / 1 \2 2 2 2 2\ 

^2\qxAq\ +<?2j ,q 2 ,m fd ,m fd ,m r ) 
Co{ql,{qi + q2) 2 1 qj,m 2 fd ,m 2 fd ,m 2 ^) = I dz 



1 l-l-Z 

dx 



x(l — x)qf + z(l — z)q\ + 2xz(qi ■ q 2 ) — (1 — z)m 2 fu — zm 2 d ' 

(A23) 

xz 



x(l — x)q\ + z(l — z)q 2 + 2xz(q\ ■ q 2 ) — (1 — z)m? L — zml ' 

(A24) 



1 pl-Z 

dz / dx 



Jo 



x(l - x)q\ + z(l - z)q\ + 2xz(qi • 92) - (1 - z)m 2 d - zm 2 u ' 



(A25) 



1 r\-z 

dx ■ 



x(l — x)qf + z(l — z)q% + 2xz(qi ■ q 2 ) — (1 — z)m d — zm fu 

Jj J j 

(A26) 



By applying the above formulae to Eqs. (|A15[) and (|A16[) , taking the first order approximation in qi, and integrating 
by x, we obtain Eq. ([T5|) . 
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